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Cn| I Abstract. For a broad range of parameter values, braneworld models display a 

remarkable property which we call cosmic mimicry. Cosmic mimicry is characterized 
by the fact that, at low redshifts, the Hubble parameter in the braneworld model is 
virtually indistinguishable from that in the LCDM (A + Cold Dark Matter) cosmology. 

■ An important point to note is that the f2 m parameters in the braneworld model and 
| in the LCDM cosmology can nevertheless be quite different. Thus, at high redshifts 

(early times) , the braneworld asymptotically expands like a matter-dominated universe 
, with the value of fi m inferred from the observations of the local matter density. At low 

redshifts (late times), the braneworld model behaves almost exactly like the LCDM 
• model but with a renormalized value of the cosmological density parameter 17^ CDM . 

The value of f^ CDM is smaller (larger) than f2 m in the braneworld model with positive 
O ' (negative) brane tension. The redshift which characterizes cosmic mimicry is related 

+3 . to the parameters in the higher-dimensional braneworld Lagrangian. Cosmic mimicry 

is a natural consequence of the scale-dependence of gravity in braneworld models. 
The change in the value of the cosmological density parameter (from £l m at high z 
• i— i , to f2^ CDM at low z) is shown to be related to the spatial dependence of the effective 

' gravitational constant G c s in braneworld theory. A subclass of mimicry models lead 

■ to an older age of the universe and also predict a redshift of reionization which is lower 
than Zrekm — 17 in the LCDM cosmology. These models might therefore provide a 
background cosmology which is in better agreement both with the observed quasar 
abundance at z > 4 and with the large optical depth to reionization measured by the 
Wilkinson Microwave Anisotropy Probe. 



PACS numbers: 04.50.+h, 98.80.Hw 
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1. Introduction 



Braneworld models have been extensively applied to cosmology, where they demonstrate 
qualitatively new and very interesting properties (see [1, 2] for recent reviews). Theories 
with the simplest generic action involving scalar-curvature terms both in the bulk and 
on the brane are not only good in modeling cosmological dark energy but, in doing 
so, they also exhibit some interesting specific features, for example, the possibility of 
superacceleration (super negative effective equation of state of dark energy w c s < — 1) 
[3, 4] and the possibility of cosmological loitering even in a spatially flat universe [5]. 

It was noted some time ago [6] that the cosmological evolution in braneworld 
theory, from the viewpoint of the Friedmann universe, proceeds with a time-dependent 
gravitational constant. In this letter, we further study this property and report another 
interesting feature of braneworld cosmology, which we call "cosmic mimicry." It turns 
out that, for a broad range of parameter values, the braneworld model behaves exactly 
as a LCDM (A + Cold Dark Matter) universe with different values of the effective 
cosmological density parameter f2 m at different epochs. This allows the model to share 
many of the attractive features of LCDM [7, 8] but with the important new property 
that the cosmological density parameter inferred from the observations of the large- 
sale structure and cosmic microwave background (CMB) and that determined from 
neoclassical cosmological tests such as observations of supernovae (SN) can potentially 
have different values. 

An important feature of this model is that, although it is very similar to LCDM 
at the present epoch, its departure from "concordance cosmology" can be significant 
at intermediate redshifts, leading to new possibilities for the universe at the end of the 
"dark ages" which may be worth exploring. 

We relate the "mimicry" properties of the braneworld cosmology with the properties 
of gravity in braneworld theories. In particular, we show that the change in the 
cosmological density parameter Q m as the universe evolves can be related to the spatial 
scale dependence of the effective gravitational constant in braneworld theory [9]. This 
can have important consequences for cosmological models based on the braneworld 
theory and calls for more extensive analysis of their cosmological history. 



2. Cosmic mimicry 

We consider the simplest generic braneworld model with an action of the form 

S = M 3 \f (ft-2A b )-2/ K]+ f (m 2 R-2a)+[ L(h ab ,<f)) . (1) 

. J bulk J brane J J brane ^ ' J brane 

Here, TZ is the scalar curvature of the metric g ab in the five-dimensional bulk, and R 
is the scalar curvature of the induced metric h ab = g ab — n a n b on the brane, where 
n a is the vector field of the inner unit normal to the brane, which is assumed to 
be a boundary of the bulk space, and the notation and conventions of [10] are used. 
The quantity K = h ab K ab is the trace of the symmetric tensor of extrinsic curvature 
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Kab — h c aV c n b of the brane. The symbol L(h ab , 0) denotes the Lagrangian density of 
the four-dimensional matter fields whose dynamics is restricted to the brane so that 
they interact only with the induced metric h ab . All integrations over the bulk and brane 
are taken with the corresponding natural volume elements. The symbols M and m 
denote the five- dimensional and four-dimensional Planck masses, respectively, Ab is the 
bulk cosmological constant, and a is the brane tension. 

Action (1) leads to the Einstein equation with cosmological constant in the bulk: 

Gab + ^b9ab = , (2) 

while the field equation on the brane is 

m 2 G ab + ah ab = r ab + M 3 (K ab - h ab K) , (3) 

where r ab is the stress-energy tensor on the brane stemming from the last term in action 
(!)• 

The cosmological evolution on the brane that follows from (2) and (3) is described 
by the main equation [3, 11, 12, 13] 



a? 3m 2 £ 2 



1 ± \ 1 + £ 2 



p + o- _K_C_ 

3m 2 6 a 4 



(4) 



where C is the integration constant describing the so-called "dark radiation" and 
corresponding to the black-hole mass of the Schwarzschild-(anti)-de Sitter solution in 
the bulk, H = a /a is the Hubble parameter on the brane, and the term n/a 2 corresponds 
to the spatial curvature on the brane. The length scale £ is defined as 

The "±" signs in (4) correspond to two different branches of the braneworld 
solutions [3]. Models with the lower ("— ") sign were called Brane 1, and models with the 
upper ("+") sign were called Brane 2 in [3], and we refer to them in this way throughout 
this paper. 

In what follows, we consider a spatially flat universe [k = 0) without dark radiation 
(C = 0). It is convenient to introduce the dimensionless cosmological parameters 

o - Po O - - O - 1 O _ A t> (c\ 

m ~3m 2 H 2) a ~3m 2 H 2 ' l ~i 2 H 2 ' Ab ~ 6H 2 ' [ ) 

where the subscript "0" refers to the current values of cosmological quantities. The 

cosmological equation with the energy density p dominated by dust-like matter can now 

be written in a transparent form: 

= n m (i+z) 3 + + 2n e ± 2^ y/n m (i+z) 3 + n a + n e + n Ah . (7) 

-"0 

The model satisfies the constraint equation 



Q m + Q a ± 2^ yjl + Q Ab = l (8) 

reducing the number of independent Q parameters. The sign choices in Eqs. (7) and (8) 
always correspond to each other if 1 + Q\ h > Q e . In the opposite case, 1 + Q\ h < Qi , 
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both signs in (8) correspond to the lower sign in (7), and the option with both upper 
signs in Eqs. (7) and (8) does not exist (see the appendix of the first paper in [3] 
for details). The signs in (8) correspond to the two possible ways of bounding the 
Schwarzschild-(anti)-de Sitter bulk space by the brane [11, 13]. 

In the "normal" case 1 + Q Ah > , substituting fl a from (8) into (7), we get 
H 2 (z) 



H 2 



(9) 



In the "special" case 1 + £l\ h < , the Brane 1 model [the lower sign in (7)] 
corresponds to both signs of (8), and the cosmological equation reads as follows: 

H\z) 



H, 



2 



n m (i+z) 3 + 1 - ft m + 2n e t 2^Jn e ^i + n Ab 



/a, 



;i+z) 3 -o m + (\ST/i + it) • (io) 



For sufficiently high redshifts, the first term on the right-hand side of (9) or (10) 
dominates, and the model reproduces the matter-dominated Friedmann universe with 
the density parameter Q m . Now we note that, for the values of z and parameters Q\ b 
and Q e which satisfy 

Q m (i+zf < ^i + n AbT ^j 2 , (ii) 

both Eqs. (9) and (10) can be well approximated as 

^ ~ n m (i+*) 3 + 1 - n m - — ^ 



~n m (i+*) 3 + i-n m - — v 7_ ^ m (i+^) 3 -^ m . (12) 



We introduce the positive parameter a by the equation 



l+^A b 



a = ^MT- (13) 



^™ = ^Tfi m , (14) 



Then, defining a new density parameter by the relation 

we get 



tf 2 (z) 



QLCDM (l + z) 3 + 1 _ Q LCDM ; (15) 



if, 



which is precisely the Hubble parameter describing a LCDM universe. [Note 
that the braneworld parameters Qe and ^A b have been effectively absorbed into a 
"renormalization" of the matter density Q m — > f2^ CDM , defined by (14). Since the value 
of the parameter f2^ CDM should be positive by virtue of the cosmological considerations, 
in the case of the upper sign in (12) and (14), we restrict ourselves to the parameter 
region a > 1, i.e., we exclude the special case of Brane 1 model with the upper sign in 
(8) from consideration.] 
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Thus, our braneworld displays the following remarkable behaviour which we refer 

to as "cosmic mimicry": 

• A Brane 1 model, which at high redshifts expands with density parameter fz m , at 
lower redshifts masquerades as a LCDM universe with a smaller value of the density 
parameter. In other words, at low redshifts, the Brane 1 universe expands as the 
LCDM model (15) with tt^ DM < fi m [where tt^ DM is determined by (14) with the 
lower ("+") sign]. 

• A Brane 2 model at low redshifts also masquerades as LCDM but with a larger value 
of the density parameter. In this case, fz m CDM > Q m with fz m CDM being determined 
by (14) with the upper ("— ") sign. 

The range of redshifts over which this cosmic mimicry occurs is given by < z <C 
z m , with z m determined by (11). Specifically, 

-\2/3 



Z m = ^ n P l/3 " 1 > ( 16 ) 

which can also be written as 

for both braneworld models. 

Examples of cosmic mimicry are shown in Fig. 1 for the Brane 1 model, and in 
Fig. 2 for the Brane 2 model. One striking consequence of Fig. 2 is that a low-density 
(fz m = 0.04) universe consisting entirely of baryons mimics a higher- density LCDM 
model (f2 m CDM = 0.2) and can therefore be in excellent agreement with the SN data. 

In view of relation (17), it is interesting to note that we can use the equations derived 
in this paper to relate the three free parameters in the braneworld model: {Qe , ^A b , ^m} 
to jfi m , z m , £l^ DM }. These relations (which turn out to be the same for Brane 1 and 
Brane 2 models) are: 

l + fi Ab _fi m CDM 



1 + Z m y , (18) 



m 



J^LCDM 



/ QLCDM O 
' m / J 'm 



^ m V fi m CDM 



vl + ^m) 3 - (19) 



Furthermore, if we assume that the value of f2 m CDM is known (say, from the analysis 
of SN data), then the two braneworld parameters Qe and can be related to the two 
parameters Q m and z m using (18), so it might be more convenient to analyze the model 
in terms of Q m and z m (instead of fig and ^A b ). 

We also note that, under condition (11), the brane tension a, determined by (8), is 
positive for Brane 1 model, and negative for Brane 2 model. 

Since the Hubble parameter in braneworld models departs from that in LCDM 
at intermediate redshifts (z > z m ), this could leave behind an important cosmological 
signature especially since several key cosmological observables depend upon the Hubble 
parameter either differentially or integrally. Examples include: 
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Q m = 1.0, z m = 46 
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0_ = 1.0, z_ = 9 



LCDM, = 0.3 



loq(1 + z) 



Figure 1. An illustration of cosmic mimicry for the Branel model. The Hubble 
parameter in three high-density Brane 1 models with 17 m = 1 is shown. Also shown 
is the Hubble parameter in the LCDM model (dotted line) which closely mimics 
this braneworld but has a lower mass density S1^ CDM = 0.3 (ft a = 0.7). The 
brane matter density (£l m ) and the matter density in LCDM are related through 



Q 



so that Q nl > f^ CDM . The redshift interval 
during which the braneworld masquerades as LCDM, so that ^branei = -^lcdm for 
z <§; z m , is z m = 9, 46, 216 (left to right) for the three braneworld models. 



the luminosity distance oIl{z): 
dj(z) f z dz' 



c 



1 + z 'Jo H(z') ' 
the angular-size distance 

c r z dz' 

d A (z) 



(20) 



(21) 



1 + z Jo H(z') ' 

the product dA_{z)H(z), which plays a key role in the Alcock-Paczynski anisotropy 
test [14], 

the product d\(z)H~ l (z), which is used in the volume-redshift test [15], 
the deceleration parameter: 
, , H'(z) 



q{z) 



H(z) 



1 + z -1 



the effective equation of state of dark energy: 

2q(z) - 1 



w(z) 



3 [1 - n m {z)) 



Hp 

H(z) 



1 + z) 



(22) 



(23) 
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a- 



LCDM, = 0.2 




0„ = 



0„ = 



0.04, 7 = 8 

0.04, z m = 32 

0.04, z m = 130 

m 



loq(1+z) 



Figure 2. An illustration of cosmic mimicry for the Brane2 model. The Hubble 
parameter in three low-density Brane 2 models with f2 m = 0.04 is shown. Also 
shown is the Hubble parameter in a LCDM model (dotted line) which mimics this 
braneworld but has a higher mass density fi|f DM = 0.2 (Q A = 0.8). The brane matter- 
density parameter (f2 m ) and the corresponding parameter of the masquerading LCDM 

ft!i CDM x 



model are related as fL 



so that f2 m < O 



LCDM 



i-y/n e /(i + n Ab ) 

The redshift interval during which the braneworld masquerades as LCDM, so that 
#brane2 = -Hlcdm f° r z *C z m , is z m — 8, 32, 130 (left to right) for the three 
braneworld models. 



the age of the universe: 

t(z) = 

the "statefinder pair" [16]: 
a 



dz' 



l+z')H(z' 



(24) 



= 1 



aH 3 
r — 1 



H" (E ,s 



(25) 



3(9-1/2)' 

the electron-scattering optical depth to a redshift z re ion 

T(Z — ) = C io {\ + z)H{zY (26) 
where n e is the electron density and &t is the Thomson cross-section describing 
scattering between electrons and CMB photons. 
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Mimicry LCDM model has fi = 0.3 



= 0.2, t = 4 

m ' m 

n m = o.i, z m = 4 

m ■ m 

fi = 0.04, z = 4 



10 



20 30 
redshift z 



40 



50 



Figure 3. The age of the universe in the Brane2 model is shown with respect to the 
LCDM value. The mimicry redshift (17) is z m = 4 so that H^ ranc (z) — #lcdm(z) 
at z <C 4. The braneworld models have Q m — 0.2, 0.1, 0.04 (bottom to top) whereas 
qlcdm = Q _ 3 _ Note that the braneworld models are older than LCDM. 



A degree of caution should be exercised when comparing the late-time LCDM 
behaviour (15) of our model with different sets of observations, since the parameter 
f2^ CDM , residing in (15), which is effectively used in determinations of the luminosity 
distance (20) and angular-size distance (21), may very well be different from the value of 
Q m inferred from observations of gravitational clustering. These issues should be kept in 
mind when performing a maximum-likelihood analysis using data belonging to different 
observational streams. 

Cosmological tests based on the luminosity distance and angular-size distance 
typically probe lower redshifts z < 2. Therefore, if the mimicry redshift is z m > 2, 
the braneworld model will, for all practical purposes, be indistinguishable from the 
LCDM cosmology on the basis of these tests alone. However, tests which probe higher 
redshifts should be able to distinguish between these models. For instance, since 
H(z) < i^LCDM(-z) in Brane2 at redshifts larger than the mimicry redshift, it follows that 
the age of the universe will be greater in this model than in the LCDM cosmology. This 
is illustrated in Fig. 3 for three distinct values of the cosmological density parameter: 
tt m = 0.2, 0.1, 0.04, all of which are lower than Q^ DM = 0.3. 

Since the late-time evolution of the universe is 

t(z) ^ —^=(1 + z)~ 3 / 2 , (27) 
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one finds, for z ^> 1, 



^LCDM 

Since f2 m < f2^ CDM in the Brane2 model, we find that the age of a Brane2 universe is 
greater than that of a LCDM universe. (The reverse is true for the Brane 1 model, for 
which Q m > ft^CDM ) 

The altered rate of expansion in the braneworld model at late times (z > z m ) also 
affects other cosmological quantities including the redshift of reionization which, for the 
Brane 2 model, becomes smaller than that in the LCDM cosmology. This is because 
the lower value of H(z) in the Brane 2 model (relative to the LCDM model), when 
substituted to (26), gives a correspondingly lower value for z Tcion for an identical value 
of the optical depth r in both models. (In fact, it is easy to see that, for the Brane 2 
model, the value of z rcion decreases with decreasing z m and Q m .) Hence, braneworld 
cosmology makes it possible to have a value for the reionization redshift which is lower 
than z ~ 17 inferred for the LCDM model from the data of the Wilkinson Microwave 
Anisotropy Probe (WMAP) [17]. 

Both an increased age of the universe and a lower redshift of reionization are 
attractive properties of the braneworld model which, as we have seen, mimics the LCDM 
cosmology at lower redshifts z < z m .§ It is important to note that the presence of high- 
redshift quasi-stellar objects (QSO's) and galaxies at redshifts z > 6 indicates that the 
process of structure formation was already in full swing at that early epoch when the 
LCDM universe was less than a billion years old. Most models of QSO's rely on a 
central supermassive black hole (Mbh <^ 1O 9 M ) to power the quasar luminosity via 
accretion. Since structure forms hierarchically in the cold dark matter scenario, the 
presence of such supermassive black holes at high redshift suggest that they formed 
through an assembly mechanism involving either accretion or mergers or both. It is not 
clear whether either of these processes is efficient enough to assemble a large number 
of high-redshift QSO's in a LCDM cosmology [18, 19]. (At the time of writing, over 
400 quasi-stellar objects with redshifts > 4 have been reported in the literature.) In 
addition, it is interesting to note the discovery of an "old" quasar APM 08279+5255 at 
redshift z — 3.91, which appears to have an older age (2.1 Gyr) than the LCDM model 
(at that redshift) and may, therefore, be problematic for concordance cosmology as well 
as several other dark-energy models [20] . Since Brane 2 can be considerably older than 
LCDM, this makes it easier to account for the existence of APM 08279+5255 within 
the braneworld context than in LCDM. (We shall return to this issue in a companion 
paper.) 

Another "surprise" for the concordance cosmology emerges from the WMAP data 
which suggest that the optical depth to reionization is r = 0.17 ± 0.06. For the LCDM 
model, this translates into a rather high redshift of reionization z = 17+5 [17], suggesting 

§ Note that the decreased redshift of reionization and the increased age of the universe are properties 
that the Brane 2 model shares with the loitering universe discussed in [5] . 



/ oLCDM 

" . (28) 
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redshift z 



Figure 4. The effective equation of state (solid line) and the deceleration parameter 
(dashed line) of the Brane2 model are shown. (The dotted line shows w = — 1 
which describes the LCDM model.) The mimicry redshift (17) is z m = 4 so that 
#branc(z) - #lcdm(z) at z < 4. The braneworld has Ct m = 0.2 whereas ^ CDM = 0.3. 



that population III stars and/or mini-QSO's were already in place at z > 17 in order 
to have efficiently reionized the universe. Whether both the structure-assembly issue 
and the high reionization redshift can be successfully accommodated within the LCDM 
paradigm remains to be seen [21]. However, both issues are important enough to compel 
the serious theorist to look for alternative models which, while preserving the manifold 
strengths and successes of the LCDM model at low z, will also be able to flexibly 
accommodate the striking properties of the universe at higher redshifts. As we have 
demonstrated in this paper, braneworld cosmology may successfully alleviate some of 
the tension currently existing between theory and observations at moderate redshifts, 
while allowing the universe to be "LCDM-like" at the present epoch. 

The effective equation of state and the deceleration parameter of the Brane 2 model 
are shown in Fig. 4. The braneworld has fl m = 0.2 and, at z < 4, masquerades as a 
higher- density LCDM model with f2]^ CDM = 0.3. Note that the effective equation of 
state (23) is a model- dependent quantity, involving the model-dependent cosmological 
parameter Q m in its definition. In our case, we use the braneworld theory as our model 
with Q m defined in (6), and the effective equation of state (23) is then redshift- dependent 
even during the mimicry period when _ff branc (f2 m , z) ~ if LCDM (f2^ CDM , z). A theorist 
who is unaware of the possibility of cosmic mimicry, when reconstructing the cosmic 
equation of state from (23) with f2^ CDM = 0.3 in the place of f2 m , will arrive at a 
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different conclusion that w — — 1. This example demonstrates some of the pitfalls 
associated with the cosmological reconstruction of the equation of state, which depends 
on the underlying theoretical model and for which an accurate knowledge of Q m is 
essential; see [16, 22, 23, 24] for a discussion of related issues. 

In concluding this section, we would like to mention an additional important 
consequence of braneworld cosmology, namely that the departure of the brane Hubble 
parameter from its LCDM counterpart is likely to affect the growth of density 
perturbations. From the linearized perturbation equation 

5 + 2H5 - AnGp 5 = 0, (29) 

where 5 = (p — p)/p is the relative density perturbation, it follows that a large value of 
H will, via the second term on the left-hand side, dampen the growth of perturbations, 
while a small value of H should hasten their growth. Although this argument is also 
likely to carry over to the braneworld, the direct application of (29) to the study of 
gravitational instability (as well as CMB anisotropy) should be treated with some 
caution since (29) was derived within the general-relativistic framework, and it is likely 
to require some modification before it can be applied to the study of cosmological 
perturbations on the brane. Indeed, according to the analysis of [4, 25], the linearized 
equation describing the evolution of inhomogeneities in the braneworld models of Dvali- 
Gabadadze-Porrati (DGP) type [26], corresponding to A b = 0, on spatial scales much 
shorter than i has the form 

5 + 2H5-AnGp(l + ^-)5 = 0, (30) 
where the time-dependent parameter /3* is given by [4] 

m = i T m[\ + ^ f )i . (si) 

Another equation which might possibly describe the behaviour of density 
perturbations in the braneworld models under consideration in this paper will be 
discussed in the next section. 

3. Properties of braneworld gravity 

3.1. Multiscale gravity 

Here, we would like to discuss the results of the previous section in the light of some 
generic properties of braneworld gravity. 

Cosmology determined by action (1) has two important scales, namely, the length 
scale (5) 

<32) 

which describes the interplay between the bulk and brane gravity, and the momentum 
scale 
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which determines the role of the brane tension in the dynamics of the brane. In a model 
characterized by the Randall-Sundrum constraint [27] 

A » s = T + 3^ = ' (34) 



the absolute value of k is equal to the inverse warping length £ warp = y — 6/Ab of the 
Randall-Sundrum solution. Note that k is negative when a < 0. We take M to be 
positive in our paper, since, in the opposite case, the massive Kaluza-Klein gravitons 
become ghosts [28, 29]. 

Following the procedure first employed in [30] for the Randall-Sundrum (RS) model 
[27] and subsequently applied in [6] to the model under consideration, we contract once 
the Gauss identity 

Rab/ = hjh b 9 h c k h d 3 n fgk j + K ac K b d - K bc K a d (35) 

on the brane and, using Eq. (2), obtain the equation 

Gab + Kfthab = 87lG c f[T ab + ^ (yJ^Qab - W a bj , (36) 



where 



3M 6 

is an important dimensionless parameter, 



SvrGW = -f^ ■ — (39) 



A eff = ^ (38) 

is the effective cosmological constant on the brane, 

I 1 

13+1 m 

is the effective gravitational constant, 

Qab = ~EE ab — E ac E c b + - (^E cd E cd — ^E 2 ^J h ab (40) 

is a quadratic expression constructed from the "bare" Einstein equation E ab = m 2 G ab — 
r ab on the brane, and E = h ab E ab . The symmetric traceless tensor W ab = n c n d W ac bd 
in (36) is a projection of the bulk Weyl tensor W a bcd- It is related to the tensor Q ab 
through the conservation equation 

D a (Q ab -M 6 W ab )=0, (41) 

where D a denotes the covariant derivative on the brane associated with the induced 

metric hab- 
it is important to note that all couplings in Eq. (36), including the effective 

cosmological and gravitational constants, are inversely proportional to (3 + 1, which 

indicates that the theory becomes singular for the special case f3 — — 1 (see [29, 31]). 
In the absence of the curvature term on the brane (m — 0), we obtain Eq. (36) in 

which 87rG e ff = 2<r/3M 6 is the gravitational constant in the RS model [27], and (3 = 0; 
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in this form, Eq. (36) was first derived in [30]. The conditions a — and Ab = are 
characteristic of the DGP model [26], which also has f3 — 0. In this model, the effective 
gravitational constant (39) turns to zero, i.e., the term linear in the stress-energy tensor 
on the brane vanishes in Eq. (36). 

Let us now turn to the "cosmic mimicry" exhibited by braneworld cosmology and 
show how it can be related to the gravitational properties of braneworld theory. In 
this context, it is remarkable that the parameter j3 introduced in (37) is very close (in 
absolute terms) to the parameter a introduced in (13) in our discussion of mimicry 
models. Specifically, 

" = i ir T0 " (42) 

so that 

\(3\ w a (43) 

when |1 — Q m \ <C Qe . This last inequality follows naturally from condition (11) for values 
of a of order unity, which are of interest to us in this paper. As a consequence, the term 
which appears in the "renormalization" of the cosmological mass density Q m in (14) is 
almost identical to the term which redefines the gravitational constant in (39). This 
coincidence can be explained by inspecting the brane equation (36). First, we note that 
cosmological solutions without dark radiation are embeddable in the anti-de Sitter bulk 
spacetime, so that W a b = for these solutions. For large cosmological matter densities, 
the quadratic expression (40) dominates in Eq. (36), and the universe is described by 
the "bare" Einstein equation m 2 G a b — r ab = 0, with the effective gravitational coupling 
being equal to 1/m 2 . As the matter density decreases, the role of this quadratic term 
becomes less and less important, and the effective gravitational coupling eventually is 
determined by the linear part of Eq. (36), i.e., by the gravitational constant (39). Thus, 
comparing (39) and (14), one has the following natural relation, valid to a high precision 
in view of (43): 

LCDM = STTGeffA, ( } 

Now we consider some relevant properties of braneworld gravity. The expression 
for Qat, in Eq. (40) is quadratic in the curvature as well as in the stress-energy tensor. 
On the other hand, the tensor W a b is related to Q a b through the conservation equation 
(41). One might, therefore, expect that the term in the parentheses on the right-hand 
side of Eq. (36), namely Q a b/M 6 — W a b, will be insignificant on sufficiently large length 
scales, and that the braneworld theory on those scales should reduce to Einstein gravity 
with the effective constants given by (38) and (39). This expectation is borne out by 
a detailed analysis [29] carried out for a positive-tension brane (er > 0) in the specific 
case when the braneworld satisfies the RS constraint (34). In this case, the gravitational 
potential of a unit mass on large scales (on the positive-tension brane) has the Newtonian 
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form with a small RS correction [29] ||: 



V(r) 



1 + 



kr > 1 , 



(45) 



3((3+l)(kr)\ 

where G e g is given by (39). 

On smaller spatial scales, kr <C 1,(3, the potential in linear theory again has the 
Newtonian form, this time with a small logarithmic correction: 



V{r) = 





kr < (46) 



and with a different expression for the effective gravitational constant [29] 



1 + 



1 



3(1 + (3) 



1 



8nm 2 



1 + 



3/3 



G, 



eft 



(47) 



0), this reproduces the result obtained for the DGP model 



For k — > (hence, also (3 
in [9] on scales r <C I. 

It is worth noting that gravity on these smaller scales kr <C l,/3, in general, 
involves the massless scalar radion, i.e., it is of scalar-tensor type. As a consequence, 
for the spherically symmetric solution, it violates the property h 00 (r) = — h~^{r), or, 
in the linear approximation, 7oo(^) = lrr{r), where 7 a /?( r ) are the components of metric 
perturbation in the spherically symmetric coordinate system. Specifically, in the model 
with the RS constraint one can obtain the relation: 

A7 _ 700 - 7rr _ 1 
700 



(4* 



700 1 + 3/3/4 ' 

Since there are stringent experimental upper bounds [32] on the left-hand side of (48) in 
the neighbourhood of the solar system (it should not exceed 10~ 5 by order of magnitude), 
if solution (46) were applicable in this domain, it would imply that only very large 
values of (3 are permissible in the braneworld theory under consideration [namely, the 
braneworld model (4) with the RS constraint (34)]. It is important to note, however, 
that the applicability region of the linear approximation (46) is bounded from below 
by a length scale which depends upon the mass of the central source, as has been 
demonstrated for the DGP model in [9]. Specifically, the dynamics of the radion develops 
strong nonlinear corrections on sufficiently small scales, leading to the breakdown of 
linearized theory. (This also creates the so-called strong-coupling problem in the DGP 
model [33].) In this case, in order to study gravity at small distances from the source, 
one must turn to the fully nonlinear theory. 

Next, let us determine the distances at which the linearized theory, leading to (46), 
breaks down and to establish the correct behaviour of the potential on such scales. In 
order to do this, we turn to the effective equation (36). Taking the trace of Eq. (36), 
we get the following closed scalar equation on the brane: 

Q 



-R + 4A 



cS 



lnG cS T 



(49) 



{(3 + 1)M 6 ' 

|| Here, we present the results for the two-brane model as the second (negative- tension) brane is taken 
to spatial infinity. 
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where the left-hand side contains terms which are linear in the curvature and in the 
stress-energy tensor while the right-hand side contains the quadratic term Q = h ab Q ab . 

Suppose that we are interested in the behaviour of gravity in the neighbourhood 
of a spherically symmetric source with density p s , total mass Ai s , and radius r s . First 
of all, we assume that one can neglect the tensor W a b and the effective cosmological 
constant in the neighbourhood of the source. As regards the effective cosmological 
constant, this assumption is natural. Concerning the tensor W a b, its smallness in the 
neighbourhood of the source represents some additional condition on the spherically 
symmetric solution. We believe that a condition of this sort is likely to arise in any 
consistent and viable braneworld theory since, without it, one has a large number of 
spherically symmetric solutions on the brane, many of them non-physical (see [34] for a 
comprehensive treatment in the framework of the RS model). 

Within the source itself, we have two qualitatively different options: an approximate 
solution can be sought either neglecting the quadratic part or linear part of Eqs. (36) 
and (49). We should choose the option that gives the smaller error of approximation in 
Eq. (49). In the first case, neglecting the quadratic part and the effective cosmological 
constant, we have 

G--8*G—»0 ^ _JL_._il_. (50) 

In the second case, we neglect the linear part, so that 

Q ab ^0 Kb~0 R + 8nG cS r ~ Ps . (51) 

(p + l)m 2 

The final expression on the right-hand side of (51) is smaller than the corresponding 
expression in (50) if 

Ps>(P + l) 2 ^ =► rg<rg~ 7;; 2 2 , (52) 
m z (p + \ym l 

where we used the relation M. s ~ p s r z s . Thus, we can expect that, in the neighbourhood 
of the source, on distances smaller than r* given by (52), the solution is determined 
mainly by the quadratic part Q ab in Eq. (36), which means that it respects the "bare" 
Einstein equation m 2 G ab = r ab to a high precision. This effect is sometimes described 
as the "gravity filter" of the DGP model [26], which screens the scalar graviton in the 
neighbourhood of the source making the gravity effectively Einsteinian. Some aspects 
of this interesting phenomenon are discussed in recent papers [35]. 

Expression (52) generalizes the length scale [9] of the DGP model, below which 
nonlinear effects become important, to the case of nonzero brane tension (nonzero 
P) and bulk cosmological constant satisfying the RS constraint (34). The observable 
gravitational constant on scales much smaller than r* will be given by 

87rG obs = \ . (53) 

For the Sun, the scale r* is estimated as 

10 16 km . . 

r,Srs ~03 + l)V» n V" < 54 > 
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which, for interesting values of (3 and , will be very large. The corresponding radius 
for the Earth is smaller only by two orders of magnitude. 

This, however, is not the full story. As argued in [4, 25], the gravitational potential 
of a spherically symmetric body on scales r* < r <C I is corrected by the cosmological 
expansion. Moreover, the critical scale r* becomes dependent on the value of the 
Hubble parameter, and can be different from (52) for values of i of the order of 
the Hubble length. Another effective gravitational constant appears in the equation 
for the growth of cosmological perturbations on the above-mentioned scales, as was 
mentioned in the previous section [see Eqs. (30) and (31)]. However, we note that, in 
this paper, interesting values of the parameter i are smaller than the Hubble length 
since we would like relation (11) to be satisfied for reasonably high redshifts, and this 
condition implies that the parameter , defined in Eq. (6), is significantly larger than 
unity. Consequently, the critical scale r* is given by Eq. (52). Note that gravity on 
scales r<r t , although close to Einstein gravity, is not exactly Einsteinian, and these 
deviations can be used to test braneworld theory on solar-system scales [4, 9, 25, 36]. 

The appearance of the distance scale r* given by (52) can also be justified by using 
an argument from cosmology. Imagine the central body to be formed of pressureless 
matter and to represent a part of the homogeneous and isotropic universe. Then 
the solution inside this body is uniquely described by the cosmological equations of 
Sec. 2. As the density parameter of the body exceeds the right-hand side of (11), 
its gravitational evolution is effectively governed by the "bare" Einstein equation 
m 2 G a b = Tab- The relevant condition describing this case is the inequality opposite 
to (11); it can be written in terms of the density p s as 

M 6 

Ps> (a T l) 2 —, (55) 

which, in view of relations (42) and (43), essentially coincides with condition (52). 

Although the preceding reasoning is applicable to both the positive-tension and 
the negative-tension brane, the current understanding of the braneworld gravitational 
physics supports only the positive-tension case (which corresponds to Brane 1 model in 
this paper), while the situation with the negative-tension brane remains unclear (at least, 
to the authors of this paper). From Eq. (36), one might expect that a negative-tension 
brane will show reasonable physical behaviour in the case \/3\ > 1 (note that (3 < for a 
negative-tension brane), in which the gravitational constant (39) is positive. However, 
direct calculation (along the lines of [29]) in the two-brane model with the RS constraint 
(34) shows that, in this case, the gravitational interaction between material bodies 
on large scales is dominated by the ghost-like radion, with the effective gravitational 
coupling 

Gradion = ~ ^ G cS , kr > 1 , (56) 

where G e g is given by the same expression (39). The radion-dominated gravity on 
these scales is formally attractive in the case G c g < 0, and is repulsive for G c q > 0. 
However, on smaller spatial scales /cr* < kr <C 1, \/3\, Newton's law similar to (46) is 
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reproduced with the gravitational constant given by (47), which is positive if \/3\ > 4/3. 
The gravity on these scales is of scalar-tensor character. On still smaller distances from 
the central source, r < r*, the theory may approach Einstein gravity with the effective 
gravitational constant (53). However, all these expectations remain to be verified by 
more refined calculations which generalize to the braneworld models of interest to us in 
this paper, including those with negative tension and without the RS constraint (34). 
In the following subsection, we propose a different approach to this problem, based on 
alternative boundary conditions in the brane-bulk system. 

3.2. The role of boundary conditions in braneworld models 

Many papers (see, e.g., the recent review [37]) have stressed that equation (36) is not 
closed on the brane in the sense that it contains the symmetric traceless tensor W a b 
whose behaviour on the brane is not determined by the dynamics of matter alone. As 
of now, this "nonclosure" of the equations on the brane constitutes one of the major 
problems confronting braneworld theory. Some additional information from the bulk is 
needed to determine the observable braneworld dynamics completely; one can say that 
some boundary conditions are to be specified to restrict the class of possible braneworld 
solutions. Usually, one tends to specify the "physical" boundary conditions somewhere 
in the bulk, by requiring the bulk to be singularity-free (see [37]) and/or by introducing 
one or more (regulatory) branes in the bulk. However, it seems to be difficult to 
implement this proposal in practice in a fully nonlinear theory, and it thus remains 
unclear whether such proposals will lead to a reasonable theory at all. A more practical 
approach to this problem, adopted in many papers (see, e.g., [38]), consists in making 
some reasonable assumptions about the tensor W a b on the brane. Here, we make a 
proposal to take this approach as a way of specifying the boundary conditions in the 
brane-bulk system. The proposal consists in specifying the boundary conditions on the 
brane, rather than in the bulk, in such a way as to close equation (36). 

Of all possible versions of boundary conditions of this kind, the simplest one appears 
to consist in setting 

W ah = (57) 

on the brane. This condition is very often imposed in practical calculations by neglecting 
the contribution of this term in various situations. Here, we do not discuss the physical 
meaning of this proposal in all detail but only consider its consequences very briefly. 

With the specification of the boundary condition in the form (57), equation (36) 
becomes a closed four-dimensional equation on the brane completely determining its 
dynamics, and solution in the bulk becomes of no crucial importance for an observer 
on the brane. Note that the specification (57) preserves all homogeneous and isotropic 
cosmological solutions with the only restriction that the dark radiation is zero, i.e., 
C = in Eq. (4). Thus we are not losing the cosmological properties of the braneworld 
models (except those for which the nonzero value of the dark radiation is crucial [5]); in 
particular, the cosmic mimicry phenomenon discussed in this paper remains intact. The 
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boundary condition under consideration affects only the inhomogeneous gravitational 
physics. 

It was shown in [39] that a spherically symmetric metric in the case (57) is unique 
and is the Schwarzschild metric in a broad region of parameters of the theory. Thus, 
in particular, there are no higher-dimensional corrections to the Newtonian potential, 
as is the case in the two-brane model [see Eqs. (45) and (46)], and the theory is not 
"spoiled" by the presence of the scalar gravitational degree of freedom. The small- 
scale gravitational physics in this case has the effective gravitational constant equal 
to 1/m 2 , as in the previous subsection. (From this viewpoint, the two-brane model 
discussed in the preceding subsection represents a braneworld model with a different 
type of boundary condition.) 

The linear equations for the cosmological scalar density perturbations in the 
braneworld model under consideration with the boundary condition (57) were, in 
fact, recently derived in [40]. Condition (57) in [40] was assumed as a reasonable 
approximation on sufficiently small spatial scales while, in this paper, we propose to 
take it as an exact boundary condition. The equation governing the evolution of the 
relativistic potential $ was shown in [40] to have the form 

^V 2 $ - 3H 2 <S> - 3HQ = -87iG* cS 5p , (58) 



where 



a' 



8nG* cS = — 2 
m 2 



( \ 

and the last equality is valid for the case of dust. As in the discussion in the end 
of the previous section, again the only change is the appearance of a time-dependent 
effective gravitational constant. Equation (58) was used to test the theory against CMB 
observations in [40]. 

We note that the proposal of the boundary conditions (57) works here for branes 
with positive as well as negative tensions and, hopefully, represents one of the ways of 
combining the attractive features of the mimicry cosmology of the Brane 2 model with 
the observable small-scale gravity. More analysis is required to test this proposal, and 
we are going to return to this issue in the future publications. 



4. Conclusion 



We have shown that braneworld cosmology, for a large region of parameter space, 
exhibits a property which can be called "cosmic mimicry." During early cosmological 
epochs, the braneworld behaves like a matter-dominated Friedmann universe with the 
usual value of the cosmological parameter fl m that would be inferred from observations 
of the local matter density. At late times, however, the universe evolves almost exactly 
like in the LCDM scenario but with a renormalized value of the cosmological density 
parameter f2j^ CDM . Specifically, a positive-tension Brane 1 model, which at high redshifts 
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expands with density parameter f2 m , at lower redshifts mimics the LCDM cosmology 
with a smaller value of the density parameter f2^ CDM < Q m . A negative-tension Brane2 
model at low redshifts also mimics LCDM but with a larger value of the density 
parameter f2^ CDM > Q m .^f The transition redshift between the early epoch and the 
late (mimicry) epoch is, in principle, a free parameter which depends upon constants 
entering the braneworld Lagrangian; see (16) or (17). 

The braneworld models discussed in this paper have interesting cosmological 
properties. For instance, in the case of Brane 1 (Brane2), the universe expands faster 
(slower) than in the LCDM scenario at redshifts greater than the mimicry redshift z m , 
whereas, for z < z m , H^, ranc (z) = H^cbm(z) in both models. The smaller value of the 
Hubble parameter at intermediate redshifts (z > few) in the case of Brane 2 leads to 
an older universe and also to a redshift of reionization which can be significantly lower 
than z ~ 17 inferred for the LCDM model from the WMAP data [17]. These features, 
together with the mimicry property, result in an interesting new cosmology which can 
successfully ameliorate the current tension between the concordance (LCDM) cosmology 
and the high-redshift universe. 

The effect of cosmic mimicry and the existence of two asymptotic density 
parameters Q m and fi^ CDM is a consequence of the time-dependence of the effective 
gravitational constant in braneworld theory [6] , which can be related to the well known 
property of the scale-dependence of the effective gravitational constant in braneworld 
models [9]. On large spatial scales, kr ^> 1, the braneworld model with positive brane 
tension (Brane 1) exhibits gravity with the renormalized effective gravitational constant 
(39), and we showed that this renormalization corresponds to the renormalization of 
the cosmological density parameter (44). The behaviour of the braneworld model with 
negative brane tension (Brane 2) on these scales is drastically different: it is dominated 
by the scalar radion which can even make the gravitational interaction repulsive; see 
Eq. (56). + The problem of this kind may not arise in models with stabilized radion [41] 
or in models with different boundary conditions for the brane-bulk system — but this 
issue has to be studied separately. 

Only two effective gravitational constants appear in the cosmology under 
consideration, given by (39) and (53), respectively, for low and high energy densities. 
However, in the local gravitational physics, there also appears the spatial distance (52) 
depending upon the mass of the central source, so that gravity in the range 

n<r<£ (60) 

of distances r from the source has a different value of the gravitational constant, given by 

(47) , and, moreover, typically has a scalar-tensor character manifested, in particular, in 

(48) . This may be important for the estimates of masses from the dynamics of clusters of 
galaxies and from gravitational lensing on these scales in the braneworld theory [9, 25]. 

The value of £1^ CDM is determined by (14) with the upper ("— ") sign for Brane 2 and the lower ("+") 
sign for Brane 1. 

+ This result was obtained for the brane with the RS constraint, but the presence of a small effective 
cosmological constant (38) is unlikely to remedy the situation. 
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On small distances from the central source, r <C r*, both positive-tension and 
negative-tension branes apparently behave similarly reproducing the Einstein gravity to 
a high precision with the gravitational constant 1/m 2 , which is the bare gravitational 
coupling in the braneworld action (1). However, this expectation is to be verified 
by refined calculations in braneworld models with arbitrary sign of brane tension and 
without the RS constraint (34). In this respect, we should note that the solution for 
a spherically symmetric source (the analog of the Schwarzschild and interior solution 
in general relativity) largely remains an open problem in braneworld theory (for recent 
progress in the DGP model, see [42]). 

Specification of boundary conditions is of crucial importance for the construction 
of a specific braneworld model since they determine the inhomogeneous gravitational 
physics on the brane. From the general viewpoint, we may regard boundary conditions 
as any conditions restricting the space of solutions of the brane-bulk system of equations 
(2), (3).* One of the possible proposals in braneworld theory that we put forward in 
this paper is to specify the boundary conditions on the visible brane, and we discussed 
one of the simplest such specifications, given by Eq. (57). In this case, the equations 
of the theory become completely closed on the brane and allow for both positive and 
negative brane tensions preserving all cosmological properties. Moreover, this proposal 
might be extended to the braneworld theory with timelike extra dimension preserving its 
attractive features [44] but without inheriting its nonphysical properties (the tachyonic 
character of the Kaluza-Klein gravitational modes). It should be stressed, however, that 
the homogeneous and isotropic braneworld cosmology is uniquely described by Eq. (4) 
and, apart from possible restriction on the dark-radiation constant C, is not sensitive 
to the specification of boundary conditions for system (2), (3). 

The cosmological model under consideration appears to safely pass the existing 
constraints on the variation of the gravitational constant from primordial abundances 
of light elements synthesized in the big-bang nucleosynthesis (BBN) and from cosmic 
microwave background (CMB) anisotropy [45]. The value of the gravitational constant 
at the BBN epoch in our model coincides with the value measured on small scales (53), 
and the effective gravitational constant (47) or (39) that might affect the large-scale 
dynamics of the universe responsible for the CMB fluctuations is within the uncertainties 
estimated in [45]. 

Cosmic mimicry in braneworld models is most efficient in the case of parameter 
a ~ 1, which, according to (13), implies that the two spatial scales, namely, the brane 
length scale given by (5) and the curvature scale of the bulk are of the same order: 
\J— 6/Ab . This coincidence of the orders of magnitude of completely 
independent scales can be regarded as some tuning of parameters, although it is 
obviously a mild one. 

* This can be compared to the view on boundary conditions in quantum cosmology, which are regarded 
there as any conditions suitably restricting the space of the wave functions of the universe. In particular, 
the boundary condition may consist in the universe having no boundary [43]. 



Cosmic Mimicry: Is LCDM a Braneworld in Disguise ? 



21 



Acknowledgments 

One of the authors (VS) thanks Bruce Bassett for an interesting conversation. The 
authors acknowledge support from the Indo-Ukrainian program of cooperation in science 
and technology sponsored by the Department of Science and Technology of India and 
Ministry of Education and Science of Ukraine. 

References 

[1] Maartens R, 2004 Living Rev. Rel. 7 7 [gr-qc/0312059] 

[2] Sahni V, Cosmological surprises from Braneworld models of dark energy, 2005 Preprint astro- 
ph/0502032 

[3] Sahni V and Shtanov Yu V, 2003 J. Cosmol. Astropart. Phys. JCAP11(2003)014 [astro- 
ph/0202346] 

Alam U and Sahni V, Supernova Constraints on Braneworld Dark Energy, 2002 Preprint astro- 
ph/0209443 

[4] Lue A and Starkman G D, 2004 Phys. Rev. D 70 101501 [astro-ph/0408246] 

[5] Sahni V and Shtanov Yu, 2005 Phys. Rev. D 71 084018 [astro-ph/0410221] 

[6] Maeda K, Mizuno S and Torii T, 2003 Phys. Rev. D 68 024033 [gr-qc/0303039] 

[7] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000), [astro-ph/9904398]. 

[8] V. Sahni, Dark Matter and Dark Energy, Lectures given at the 2nd Aegean Summer School on the 

Early Universe, Ermoupoli, Island of Syros, Greece, astro-ph/0403324. 
[9] Gruzinov A, 2005 New Astron. 10 311 [astro-ph/0112246] 

Porrati M, 2002 Phys. Lett. B 534 209 [hcp-th/0203014] 

Dvali G, Gruzinov A and Zaldarriaga M, 2003 Phys. Rev. D 68 024012 [hcp-ph/0212069] 
[10] Wald R M, 1984 General Relativity, (Chicago: The University of Chicago Press) 
[11] Collins H and Holdom B, Phys. Rev. D 62 105009 [hep-ph/0003173] 
[12] Shtanov Yu V, On brane-world cosmology, 2000 Preprint hep-th/0005193 
[13] Deffayet C, 2001 Phys. Lett. B 502 199 [hep-th/0010186] 
[14] Alcock C and Paczynski B, 1979 Nature 281 358 

Kujat J, Linn A M, Scherrer R J and Weinberg D H, 2001 ApJ S72 1 [astro-ph/0112221] 
[15] Newman J A and Davis M, 2000 ApJ 534 Lll (2000) [astro-ph/0109131] 

[16] Sahni V, Saini T D, Starobinsky A A and Alam U, 2003 JETP Lett. 77 201 [astro-ph/0201498] 
Alam U, Sahni V, Saini T D and Starobinsky A A, 2003 MNRAS 344 1057 [astro-ph/0303009] 

[17] Spcrgcl D N et al, 2003 Astrophys. J. Suppl. 148 175 (2003) [astro-ph/0302209] 
Tegmark M et al, 2004 Phys. Rev. D 69 103501 (2004) [astro-ph/0310723] 

[18] Richards G T et al, 2004 Astron. J. 127 1305 [astro-ph/0309274] 

[19] Haiman Z and Quataert E, The Formation and Evolution of the First Massive Black Holes, 2004 

Preprint astro-ph/0403225 
[20] Hasinger G, Schartel N and Komossa S 2002 Astrophys. J. 573 L77 

Friaga ACS, Alcaniz J S and Lima J A S 2005 Mon. Not. Roy. Astron. Soc. 362 1295 [astro- 

ph/0504031] 

Jain D and Dev A 2005 Age of High Rcdshift Objects - a Litmus Test for the Dark Energy 
Models Preprint astro-ph/0509212 
[21] Barkana R and Loeb A, In the Beginning: The First Sources of Light and the Reionization of the 
universe, 2000 Preprint astro-ph/0010468 

Ciardi B and Ferrara A, The First Cosmic Structures and their Effects, 2004 Preprint astro- 
ph/0409018 

[22] Maor I et al, 2002 Phys. Rev. D 65 123003 [astro-ph/0112526] 

[23] Alam U, Sahni V, Saini T D and Starobinsky A A, 2004 MNRAS 354 275 [astro-ph/0311364] 



Cosmic Mimicry: Is LCDM a Braneworld in Disguise ? 



22 



Alam U, Sahni V and Starobinsky A A, 2004 J. Cosmol. Astropart. Phys. JCAP06 (2004)008 
[astro-ph/0403687] 

[24] Bassett B S, Corasaniti P S and Copeland E J, 2004 ApJ 617 LI [astro-ph/0407364] 
[25] Lue A and Starkman G, 2003 Phys. Rev. D 67 064002 [astro-ph/0212083] 

Lue A, Scoccimarro R and Starkman G, 2004 Phys. Rev. D 69 044005 [astro-ph/0307034] 

Lue A, Scoccimarro R and Starkman G D, 2004 Phys. Rev. D 69 124015 [astro-ph/0401515] 
[26] Dvali G, Gabadadze G and Porrati M, 2000 Phys. Lett. B 485 208 [hcp-th/0005016] 

Dvali G and Gabadadze G, 2001 Phys. Rev. D 63 065007 [hep-th/0008054] 
[27] Randall L and Sundrum R, 1999 Phys. Rev. Lett. 83 3370 [hep-ph/9905221] 

Randall L and Sundrum R, 1999 Phys. Rev. Lett. 83 4690 [hcp-th/9906064] 
[28] Padilla A, 2004 Class. Quantum Grav. 21 2899 [hep-th/0402079] 
[29] Shtanov Yu and Viznyuk A, 2005 Class. Quantum Grav. 22 987 [hep-th/0312261] 
[30] Shiromizu T, Maeda K and Sasaki M, 2001 Phys. Rev. D 62 024012 [hcp-th/9910076] 
[31] Smolyakov M N, 2004 Nucl. Phys. B 695 301 [hcp-th/0403034] 
[32] Bertotti B, less L and Tortora P, 2003 Nature 425 374 

Will C M, Was Einstein Right? Testing Relativity at the Centenary, 2005 Preprint gr- 

qc/0504086 

[33] Luty M A, Porrati M and Rattazzi R, 2003 J. High Energy Phys. JHEP09(2003)029 [hep- 
th/0303116] 

Rubakov V A, Strong coupling in brane-induced gravity in five dimensions, 2003 Preprint hep- 
th/0303125 

[34] Visscr M and Wiltshire D L, 2003 Phys. Rev. D 67 104004 [hep-th/02 12333] 

[35] Kalopcr N, 2005 Phys. Rev. Lett. 94 181601; Erratum ibid. 95 059901 [hcp-th/0501028] 

Kaloper N, 2005 Phys. Rev. D 71 086003 [hep-th/0502035] 
[36] Iorio L, On the effects of the Dvali-Gabadadze-Porrati braneworld gravity on the orbital motion of 

a test particle, 2005 Preprint gr-qc/0504053 
[37] Lue A, The Phenomenology of Dvali-Gabadadze-Porrati Cosmologies, 2005 Preprint astro- 

ph/0510068 

[38] Maartcns R, 2003 Progr. Theor. Phys. Suppl. 148 213 [gr-qc/0304089] 

[39] Kofinas G, Papantonopoulos E and Pappa I, 2002 Phys. Rev. D 66 104014 [hcp-th/0112019] 

Kofinas G, Papantonopoulos E and Zamarias V, 2002, Phys. Rev. D 66 104028 [hcp-th/0208207] 

[40] Sawicki I and Carroll S M, Cosmological Structure Evolution and CMB Anisotropics in DGP 
Braneworlds , 2005 Preprint astro-ph/0510364 

[41] Goldberger W D and Wise M B, 1999 Phys. Rev. Lett. 83 4922 [hcp-ph/9907447] 

[42] Gabadadze G and Iglesias A Schwarzschild solution in brane induced gravity, 2004 Preprint hep- 
th/0407049 

[43] Hartlc J B and Hawking S W, 1983 Phys. Rev. D 28 2960 

[44] Shtanov Yu and Sahni V, 2003 Phys. Lett. B 557 1 [gr-qc/0208047]. 

[45] Umczu K, Ichiki K and Yahiro M, 2005 Phys. Rev. D 72 044010 [astro-ph/0503578] 



